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ABSTRACT 


The  increasing  interest  in  long  term  (inter-annual) 
weather  changes  and  their  relation  to  processes  in  the  ocean 
is  beginning  to  illuminate  the  need  for  and  the  lack  of  long 
term  records  of  physically  significant  variables  which  occur 
over  the  vast  oceanic  regions  of  the  northern  hemisphere. 
An  attempt  is  made  here  to  evaluate  the  accuracy  of  a  hind 
cast  time  series  of  surface  wind  vector  components  and  speeds 
for  the  period  of  January  1960  through  December  1969. 

The  quasi -geostrophi c  model  used  to  calculate  these 
records  is  described  as  well  as  the  twelve-hourly  surface 
pressure  data  which  were  used  as  input.   The  central  moments 
of  the  probability  distributions  of  the  computed  records  are 
compared  to  those  of  corresponding  time  series  observed  on 
Ocean  Station  Vessels.   Linear  correlation  coefficients 
between  observed  and  computed  records  were  found  to  average 
0.81  for  the  components  of  the  wind  vector  and  0.65  for  the 
wind  speed.   Regression  relations  between  computed  and 
observed  records  also  are  presented.   Spectral  analysis  of 
low  pass  filtered  records  showed  coherence  between  observed 
and  computed  records  increased  with  decreasing  frequency. 
The  accuracy  of  the  computed  records  in  low  latitudes  also  is 
estimated . 
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I   INTRODUCTION 

A.    GENERAL 

One  of  the  most  important  variables  encountered  in  the 
study  of  physical  oceanography  and  air-sea  interaction  is 
the  wind  field  at  the  sea  surface.   To  a  large  degree,  the 
surface  wind  velocity  and  its  variability,  both  in  time  and 
space,  determine  the  exchange  of  heat  and  momentum  across 
the  air-sea  interface.   Through  such  exchange  the  wind  affects 
not  only  oceanic  surface  phenomena,  such  as  waves  and  surface 
currents,  but  also  features  within  the  sea  such  as  the 
vertical  thermal  structure  and  deep  currents. 

The  smallest  scale  of  concern  in  this  study  is  the 
synoptic  scale.   Such  a  scale  covers  a  period  of  time  from 
several  hours  to  several  days  and  a  distance  from  several 
tens  of  kilometers  to  several  hundred  kilometers.   Past 
attempts  at  estimating  the  surface  wind  field  on  this  scale 
have  been  of  two  types.   One  could  either  average  wind  obser- 
vations collected  over  a  given  area  for  a  specified  period 
of  time  (Seckel,  1970)  or  one  could  estimate  the  surface  winds 
from  a  similar  collection  of  observations  of  other  variables, 
such  as  surface  atmospheric  pressure  (Roden,  1974). 

The  accuracy  of  synoptic  scale  surface  wind  fields 
derived  directly  from  wind  observations  suffers  from  two 
major  drawbacks:  a  paucity  of  observations,  and  the  high 
variability  of  the  small  scale  wind  field.   The  small  number 
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of  observations  usually  available  in  a  typical  mid-ocean 
region  over  a  three-or-four-day  period  are  very  often  grouped 
together  in  a  small  section  of  the  region,  often  along  major 
shipping  lanes  or  along  the  track  of  a  single  ship.   In 
addition,  major  portions  of  the  open  ocean  lack  any  reports 
of  wind  conditions  for  weeks  at  a  time.   As  a  result,  the 
few  reports  available  must  necessarily  be  used  to  describe 
large  areas  when  estimating  the  synoptic  scale  wind  field. 
Verploegh  (1967)  found,  however,  that  simultaneous  wind  speed 
observations  from  ships  less  than  sixty  miles  apart  in 
"synopti cal ly  homogeneous  areas"  differed  significantly.   His 
analysis  indicated  the  correlation  coefficient  between  such 
observations  was  less  than  0.7. 

Because  of  such  considerations,  the  second,  indirect 
approach  is  generally  used  to  estimate  the  synoptic  scale 
surface  wind  field.   This  approach  depends  on  the  fact  that 
the  distribution  of  surface  atmospheric  pressure  can  be  more 
accurately  described  than  can  the  surface  wind  field  from  a 
limited  number  of  observations,  due  to  the  lower  variability 
in  time  and  space  of  the  surface  atmospheric  pressure  field. 
An  accurate  estimate  of  the  synoptic-scale  surface  wind  field 
may  be  made  by  applying  the  geostrophic  approximation  to  the 
surface  atmospheric  pressure  field. 
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B.    OBJECTIVES 

Although  the  relation  of  the  geostrophic  to  observed 
wind  is  a  popular  topic,  this  study  is  not  intended  to  be  a 
detailed  examination  of  that  relation.   Rather,  it  is  limited 
to  an  evaluation  of  one  particular  surface  wind  model  used, 
in  this  study, to  hindcast  time  series  of  surface  winds  over 
oceanic  areas.   The  model  is  quasi -geostrophi c ,  incorporating 
the  effects  of  friction.   It  is  one  part  of  operational 
oceanographi c  models  used  at  Fleet  Numerical  Weather  Central 
and  at  Fleet  Weather  Central,  Rota  (Spain). 

As  a  result  of  this  study,  the  ten-year  average  and 
linear  trend  of  observed  wind  conditions  at  several  points 
in  the  North  Atlantic  and  North  Pacific  Oceans  were  determined 
and  compared  to  those  calculated  from  the  model.   The  primary 
objective  of  the  study  was  however,  to  determine  the  accuracy 
of  the  computed  time  series  and  to  determine  if  calibration 
factors  could  be  applied  directly  to  them  in  order  to  increase 
their  accuracy. 


The  model  was  developed  by  the  author  while  at  Fleet 
Numerical  Weather  Central;  most  of  the  work  of  this  study  was 
carried  out  while  the  author  was  at  the  Environmental  Prediction 
Research  Facility. 
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II.   APPROACH 

A.    THEORETICAL  BACKGROUND 

1 .    The  Geostrophi c  Wind  Equation 

Geostrophic  wind  has  been  considered  to  be  an 
acceptable  estimate  of  the  true  surface  wind  in  many  studies 
of  oceanography  and  air-sea  interaction  (Roden  p_£.  cvt .  , 
Namias,  1963).   Geostrophy  assumes  a  balance  between  the 
pressure  gradient  force  and  the  Coriolis  force  and  ignores 
acceleration,  friction,  and  vertical  motion.   This  balance  is 
gi  ven  in  equation  CO: 


Vn  =  k  X  -If  vHP 
g       pf   H 


(1) 


where  V   =  geostrophic  wind  vector 
P   =  surface  pressure 
v u  =  —  i  +  —  j   (horizontal  gradient  operator) 

n      o  X        ay 

X  =  vector  product  operator 

f  =  2oj  sin  0  (the  Coriolis  parameter) 

u  =  angular  rotation  of  the  earth 

0  =  latitude 

p  =  ai  r  densi  ty 

k  =  unit  vertical  vector 
Although  this  approximation  is  widely  used,  its  accuracy  is 
limited  by  the  required  assumptions  indicated  above. 


15 


The  true  surface  wind  vector  is  the  resultant  of 
the  geostrophic  wind  vector  and  an  ageostrophic  wind  vector. 
The  direction  and  magnitude  of  the  ageostrophic  wind  vector 
determines  the  accuracy  of  the  geostrophic  approximation. 
Unfortunately  the  ageostrophic  vector  has  been  shown  to  vary 
significantly  over  the  Northern  Hemisphere  (Roll,  1965; 
Carstensen,  1967;  Brummer,  e_t  aj_.  ,  1974).   Such  variability 
makes  the  geostrophic  wind  hard  to  relate  to  the  surface  wind 
in  studies  on  synoptic  or  larger  scales. 

The  accuracy  of  the  geostrophic  approximation  is 
also  limited,  in  practical  application,  by  the  temporal 
variability  of  the  surface  pressure  field.   In  practice,  the 
geostrophic  approximation  is  often  applied  to  surface  pressure 
distributions  averaged  over  a  month  or  more  (Namias,  op  cj_t .  ) . 
Such  a  procedure  results  in  calculating  the  vector  resultant 
of  the  geostrophic  wind  for  the  averaging  period.   Obviously, 
non-conservative  phenomena  which  depend  on  the  surface  wind 
cannot  be  treated  using  such  a  resultant  wind  field.   In  the 
study  of  such  phenomena,  the  resultant  vector  is  assumed  to 
represent  the  mean  or  steady-state  wind  prevailing  over  the 
averaging  interval. 

Whether  or  not  the  resultant  wind  is  a  good  esti- 
mate of  the  mean  wind  depends  on  the  variance  of  the  wind 
vector  over  that  interval.   The  mean  wind  speed,  WM ,  and  the 
resultant  wind  speed,  WR  are    given  below  in  terms  of  components 
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w 


T 

E 

t=l 


(u2t  ♦  v2t): 


„  -   ^ — ~  (2) 


(3) 


where  U.,  V.  are  orthogonal  components  of  the  wind  vector  at 
time  t  for  1  <_   t  <_  T .   (Averaging  twelve  hourly  data  over 

thirty  days,  we  have  1  <_  t  <_   60).   The  difference  between 

2      2 

WM  and  WR  can  be  expressed  in  terms  of  the  variance  of  each 

2      2  2 

of  the  components,  a      and  a    ,  and  of  the  wind  speed,  a    ; 


u2   „2    2,2    2 

WM  "  WR  =  %  +  av  "  aw 


(4) 


where 


T   II*   (    1      uA2 

t=i    '      yt=i    '  y 


(5) 


2        2 

and  similarly  for  a    ,    and  a  .   As  will  be  shown  in  Section 

III-A,  the  difference  expressed  in  equation  (4)  results  in 
a  drastic  under-estimate  of  wind  speed  when  the  resultant 
wind  is  used  in  place  of  the  mean  wind  speed  in  problems 
involving  the  square  of  the  wind  speed. 

2 .    The  Surface  Wind-Geostrophic  Mind  Relation 

Although  this  is  not  a  study  of  the  general  relation 
of  the  geostrophic  wind  to  the  surface  wind,  a  few  words 
regarding  this  relation  are  appropriate  for  background.   As 
indicated,  the  surface  wind  vector  can  be  resolved  into 
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geostrophic  and  ageostrophi c  components.   The  ageostrophic 
component  can  be  assumed  to  consist  of  a  frictional  component 
and  an  acceleration  component  (Briimmer,  e_t  aj_.  ,  op  ci  t .  ) . 

The  frictional  component  is  generally  accepted  as 
the  dominant  ageostrophic  component  (Verploegh,  op  c  i  t . , 
Brummer,  et  a]_.  ,  op  c  i  t .  ) .   Haltiner  and  Martin  (1957)  suggest 
that  the  frictional  component  is  dependent  on  the  stability 
of  the  surface  layer  of  the  atmosphere.   Hasse  (1974)  has 
demonstrated  that  the  dependence  of  the  surface  wind  speed  on 
the  surface  pressure  gradient  is  an  order  of  magnitude  greater 
than  its  dependence  on  stability. 

The  effects  of  stability  have  not  been  explicitly 
incorporated  into  the  model;  to  some  degree  however,  they 
have  been  implicitly  incorporated  through  the  latitude 
dependence  of  terms  in  the  various  equations  of  the  model. 
The  major  effects  of  stability  are  assumed  to  occur  in  the 
region  between  45°N  and  25°N.   This  implicit  incorporation  of 
stability  is  shown  in  the  ratio  of  the  surface  wind  speed  to 
the  geostrophic  wind  speed  in  Figure  1. 

B.    METHOD  OF  COMPUTING  THE  SURFACE  WIND 

The  calculations  were  performed  on  the  Fleet  Numerical 
Weather  Central  (FNWC)  63x63  polar  stereograph ic  grid  of  the 
northern  hemisphere  (FNWC,  1974).   On  this  square  grid  the 
equator  is  an  inscribed  circle,  i.e.,  tangent  to  the  grid 
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boundary.   The  grid  mesh  length  at  60°N  is  381  km  and  increases 
with  increasing  latitude.   The  data,  although  calculated  on  a 

polar  stereographic  grid,  will  be  discussed  relative  to  a 

2 
Mercator  grid. 

The  components  of  the  horizontal  pressure  gradient  at  a 

given  grid  point  (i,j),  were  approximated  to  the  second  order 

by  equation  (6): 


where   K 


1 


||(i  j)  =  MLLLil  Kl{P(i-2,j)  -  K2[p(i-i,j) 


-  P(i+l,j)]  -  P(1+2,j)} 

1/12 
8 


(6) 


MF(i,j) 


0 


=  381  km  (mesh  length) 

1+SIN(60°)       ,   .     ,    ..     . ,   . 
=  i +c t m  )  a  \ —  =  maP  factor  of  the  grid  at 

point  (i  ,j) 
=  latitude  of  grid  point  ( i  ,  j )  , 


g  p 
and  similarly  for  ry(i,j).   Using  (1)  with  a  constant  air 

-  3      -  3 
density  of  1.22  x  10    gm  cm   ,  the  geostrophic  wind  vector 


Both  Mercator  and  polar  stereographic  projections  are 
conformal  projections  of  the  earth  and  as  a  result,  angles 
are  preserved  when  transformed  from  one  to  the  other  (Taylor, 
1955). 
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at  each  point  was  found.    The  magnitude  of  this  vector  (the 
geostrophic  wind  speed)  was  then  reduced  as  a  function  of  the 
latitude  of  the  grid  point  according  to  (7): 


Vs(i,j)  =  A  [B(0)]  Vg(i,j) 


(7) 


where  V  (i,j)  =  surface  wind  speed  at  grid  point  (i,j) 

V  ( i  ,  j )  =  geostrophic  wind  speed  at  grid  point  (i,j) 
A        =  0.93  4 

(0.65+0.2(0/25°) ;  0  <  25° 
B(0)     =<0.85  ;  25°  <  0  <  45° 

(o. 75  +  0.1  ^9^;^;  45°  <  0 


The  resulting  ratio  of  (V  /V  ),  which  is  strongly 
dependent  on  latitude,  was  suggested  by  data  from  Roll 
(op  ci  t .  )  and  Carstensen  (op_  c  i  t .  )  .   The  value  of  (V  /V  ) 
used  in  the  model  is  plotted  in  Figure  1  along  with  values 


The  sine  function,  SILA,  used  in  the  calculation  of  the 
Coriolis  parameter  for  the  geostrophic  wind  calculation  and 
in  the  calculation  of  the  deflection  angle  (equations  1  and 
8)  was  modified  below  35°N.   This  was  done  in  an  attempt  to 
reduce  the  impact  of  this  parameter  in  the  calculation  of 
equatorial  winds.   The  modification  was  suggested  by  the 
subroutine  SNLTSRX  in  the  Fleet  Numerical  Weather  Central 
program  1 i brary . ( FNWC ,  op  c  i  t .  )  .   The  value  of  SILA  south  of 
3  5° N  is  given  by  the  tangent  to  the  sine  curve  at  3  5  °  N  ,  i.e., 
SILA(0)  -  0.0144(SIN-1  (0)  )  +  0.075  for  0  <_35°. 

4 
The  factor  A  resulted  from  erroneously  using  408  km 

for  the  grid  mesh  length  in  the  model  (D.  Dale,  FNWC, 

personal  communication). 
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found  by  Roll  and  by  Carstensen.   Although  Carstensen  derived 

his  relation  from  a  much  smaller  data  set  than  that  in  Roll's 

study,  the  major  features  of  the  latitude  dependence  are 

similar.   The  value  of  (V  /V  )  used  in  the  present  model 

s   9 

attempts  to  reproduce  the  first  order  dependence  on  latitude: 
i.e.,  the  maximum  in  the  mid-latitudes,  sharply  decreasing 
with  decreasing  latitude  toward  the  equator  and  a  more  modest 
decrease  with  increasing  latitude  above  45°N. 

After  being  reduced  in  speed,  the  surface  wind  vector 
is  then  rotated  to  the  left  of  the  geostrophic  wind  vector 
(toward  low  pressure).   Together,  deflection  and  reduction 
of  the  geostrophic  wind  vector  simulate  the  effects  of  friction 
(Haltiner  and  Martin,  op  c  i  t .  )  .   The  deflection  angle,  a, was 
cal cul ated  from  (8) : 


a  =  K1  (K2  -  K3V^)  /  (1+SIN  0) 


(8) 


where  K-,  -  1  .475 
K2  =  22.5 
K3  =  1 .75  x  10"2 
V   =  surface  wind  speed,  from  (7)  . 

Figure  2  is  a  plot  of  a  against  latitude  for  several  values 
of  surface  wind  speed. 

The  dependence  of  a  on  the  square  of  the  surface  wind 
speed  strongly  decreases  the  deflection  angle  associated  with 
high  surface  wind  speeds.   This  approach  was  based  on  the 
assumption  that  high  winds  (associated  with  well-defined 
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cyclones)  generally  show  less  cross-isobar  flow  than  do  the 
low  winds  associated  with  nearly  flat  surface  pressure 
distributions.   The  assumed  reduction  in  the  deflection  angle 
as  wind  speed  increases  could  be  caused  by  a  slight  reduction 
in  the  frictional  drag  experienced  at  higher  values  of  surface 
wind  speed.   In  this  model,  however,  this  reduction  becomes 
significant  only  at  relatively  high  wind  speed.   For  example, 
at  20  m  sec"   the  deflection  angle  is  75%  of  its  value  at 
1  rn  sec"  .   At  30  m  sec"   the  deflection  angle  has  been 


reduced  to  30%  of  the  value  at  1  m  sec 


-  i 


As  with  the  value  of  (V  /V  ),  the  latitude  dependence 
of  m  also  was  suggested  by  Roll  and  Carstensen.   Figure  3  is 
a  plot  of  the  value  of  a  against  latitude  as  found  by  Roll, 
by  Carstensen,  and  as  used  in  this  model.   The  values  of  a 
froE  the  model   are  those  for  surface  wind  speeds  of  5  and 
15  ■  sec"  .   Less  than  15%  of  the  reported  wind  speeds  in 
this  study  were  outside  this  range.   As  in  the  case  of  the 
reduction  factor,  the  expression  for  a  attempted  to  reproduce 
the  najor  features  of  the  latitude  dependence  shown  by  both 
Roll's  and  Carstensen's  data. 


C.    DATA  GENERATION 

1 .    Input  Data  for  the  Model 

The  model  was  used  to  calculate  a  ten-year  time 
series  of  the  surface  winds  at  twelve  hour  intervals  at  grid 
poiats  near  international  Ocean  Stations.   The  input  for  the 
model  consisted  of  surface  pressure  fields  on  the  FNWC  grid. 
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All  the  available  surface  pressure  analyses  for  synoptic 
periods  of  0000  GMT  and  1200  GMT  were  used  for  the  ten-year 
period  from  0000  GMT  1  January  1960  to  1200  GMT  31  December 
1969.   However,  not  all  the  analyses  for  this  period  were 
available  and  this  required  some  conditioning  of  the  output 
data  which  will  be  discussed  in  Section  II-D. 

The  number  of  surface  pressure  analyses  used  for 
each  month  in  the  ten-year  period  are  shown  in  Table  I.   In 
general,  the  analyses  were  fairly  evenly  distributed  over 
the  twelve  months.   The  month  of  April  had  the  poorest 
representation  with  488  (81%)  of  the  600  possible  analyses 
used,  while  November  was  the  best  represented  with  534  (89%) 
of  the  600  possible  analyses.   Most  of  the  missing  data  were 
concentrated  in  the  period  from  April  1960  through  June  1962. 
Within  this  period  only  one  analysis  per  day  was  available. 
As  a  result,  only  365  (50%)  of  the  730  possible  analyses  were 
used  for  1961.   The  best-represented  year  was  1967  for  which 
728  (99.7%)  of  the  730  possible  analyses  were  available.   The 
missing  analyses  totaled  1098  (15%)  of  the  7306  possible  in 
the  ten-year  interval. 

According  to  Bakun  (1973),  who  used  the  same  set 
of  surface  pressure  analyses  to  calculate  coastal  upwelling 
indices,  the  analyses  between  January  1960  and  July  1962 
were  from  the  National  Climatic  Center  at  Asheville,  North 
Carolina.   The  analyses  from  July  1962  through  December  1969 
were  produced  by  objective  computer  analyses  at  FNHC. 
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2 .  The  Observational  Data 

The  observed  wind  data  were  taken  at  six  interna- 
tional Ocean  Stations  (OS).   These  were  the  only  continuous 
time  series  of  observed  wind  speed  and  direction  available 
for  the  ten  years  of  this  study.   As  these  data  were  assumed 
to  be  of  high  quality,  no  attempt  was  made  to  determine  their 
accuracy  other  than  gross  error  checking  as  described  in 
Section  1 1  - D .   The  observations  were  thus  used  as  the  standard 
against  which  the  calculated  wind  records  were  verified. 

The  observed  wind  speed  and  direction  were  extracted 
from  the  Surface  Marine  Observation  Tape  Family  11  from  the 
National  Weather  Record  Center  for  the  following  Ocean 
Stations:   BRAVO  (OSB),  DELTA  (OSD),  JULIET  (OSJ),  PAPA  (OSP), 
VICTOR  (OSV),  and  NOVEMBER  (OSN).5   Table  II  indicates  the 
position  of  each  OS,  the  position  of  the  computational  grid 
points,  and  the  distance  between  each  OS  and  its  correspond- 
ing grid  point. 

3.  Data  Used  in  Analysis 

Both  the  observed  and  computed  data  were  transformed 
to  zonal  (East-West)  and  meridional  (North-South)  components. 
The  resulting  set  of  three  time  series  for  each  OS  and  each 
grid  point  (a  total  of  39)  consisted  of  the  values  of  zonal 


All  the  Ocean  Station  records  with  the  exception  of 
NOVEMBER  (OSN)  were  compared  to  computational  time  series 
from  a  single  grid  point  near  their  location.   Time  series 
were  calculated  at  two  grid  points  bracketing  OSN  in  order 
to  estimate  the  variability  of  the  computed  wind  over  a 
single  grid  mesh  length.   These  computation  grid  points  will 
be  referred  to,  henceforth,  as  NIC  and  N2C. 
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and  meridional  components  and  the  wind  speed  at  0000  GMT  and 
1200  GMT  for  the  ten-year  period  from  1  January  1960  to  31 
December  1969.   It  should  be  noted  here  that  a  negative  zonal 
value  indicates  a  wind  from  the  West  and  a  negative  meridional 
value  indicates  a  wind  from  the  South. 

D.    DATA  CONDITIONING 

1 .    Gaps  and  Gross  Errors 

A  check  for  gross  errors  was  performed  following 
extraction  and  transformation  of  the  observed  and  computed 
wind  vectors.   A  wind  speed  of  50  m  sec"   was  selected  as  the 
maximum  allowable  value;  this  established  a  gross  error 
criterion  in  order  to  eliminate  only  those  errors  caused  by 
hardware  failure  such  as  tape  parity.   No  wind  speed  values 
in  excess  of  50  m  sec"   were  found  in  either  the  computed  or 
observed  data. 

After  gross  error  checking  was  completed,  ensemble 
analysis  was  performed  on  the  observed  and  computed  time 
series.   Following  that,  the  gaps  in  the  records  were  removed. 
The  computed  records  differed  markedly  frora  the  observational 
data  in  both  the  number  and  length  of  gaps.   For  the  eighteen 
observational  records,  the  number  of  single  point  gaps  per 
record  averaged  about  200,  while  about  300  points  per  record 
were  missing  in  gaps  of  two  or  more  consecutive  points.   Only 
about  150  points  per  record  were  missing  from  the  computed 
time  series  in  the  form  of  multiple  point  gaps.   However,  for 
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these  computed  records  there  were  about  950  single  point 
gaps,  due  primarily  to  missing  pressure  analyses  within  the 
period  of  April  1960  through  June  1962. 

Because  of  the  high  proportion  of  single  missing 
data  points  in  the  calculated  record,  a  single  missing  value 
was  replaced  by  linear  interpolation  between  the  preceding 
and  following  value.   For  gaps  consisting  of  two  or  more  data 
points,  the  missing  data  were  replaced  by  the  time  series 
average.   This  procedure  was  suggested  as  a  means  of  filling 
gaps  while  not  introducing  spurious  data  into  the  records 
which  would  be  subjected  to  subsequent  spectral  analysis. 
2.    Filtering 

Filtering  was  done  on  the  computed  and  observed 
records  in  order  to  allow  comparison  of  records  with  nearly 
equal  data  content.   Because  of  the  substantial  missing  data 
in  the  high  frequency  part  of  the  computed  record  (due  to 
single  point  gaps)  and  because  of  the  stated  interest  in 
large  scale  phenomena,  a  low-pass  filter  was  applied  to  each 
of  the  computed  and  observed  time  series  after  the  gaps  were 


E.  Thornton,  Dept.  of  Oceanography,  Naval  Postgraduate 
School,  personal  communication. 
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removed.   A  five-point  binomial  filter  indicated  in  equation 
(9)was  used. 


Y(t)  =  s   a.  •  X(t+lAt) 
i=-2   n 


(9) 


where   Y(t) 


ai 


X(t+iAt) 

At 


=  filtered  value  at  time  t 

=  binomial  coefficient  (0.38,  0.25,  0.06, 
for  i  =  0,  + 1 ,  +2) 

=  unfiltered  value  of  the  time  series  at 
time  t+iAt 

=  sample  interval  (0.5  days) 


Figure  4  shows  the  effect  of  applying  this  filter  to  a  typical 
segment  of  one  of  the  computed  time  series.   The  amplitude  of 
the  residual  series  (unfiltered  minus  filtered)  is  approxi- 
mately the  same  amplitude  as  the  filtered  series,  indicating 
that  approximately  50%  of  the  variance  of  the  unfiltered 
series  was  removed  by  filtering.   This  estimate  of  50%  was 
confirmed  by  a  comparison  of  the  filtered  and  unfiltered  record 
variance  (not  shown).   Before  any  filtering  or  gap  removal  was 
done,  however,  the  observed  and  computed  records  were  analyzed 
as  ensembles  in  order  to  determine  the  relation  between  them. 


Such  a  filter  has  a  frequency  response,  R(f),  at 
frequency  f  , 

R(f)  -  cos4  (irfAt)  (10) 

(Panofsky  and  Brier,  1958). 
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III.   RESULTS 

A.    ENSEMBLE  ANALYSIS 

In  this  section  the  ensembles  of  the  observed  and 
computed  data  are  compared.   The  term  ensemble  is  taken  to 
mean  the  set  of  values  which  constitute  a  given  time  series 
of  the  variable  in  question.   Due  to  the  length  of  the  records 
(ten  years),  significant  questions  about  short-term  accuracy 
of  the  model  cannot  be  treated  by  ensemble  analysis.   The 
problem  of  short-term  accuracy  is  investigated  through 
spectral  analysis  as  reported  in  Section  III-B. 

1 .    Central  Moments 

Both  the  observed  and  computed  data  were  considered 

to  have  come  from  specific  probability  densities.   While  these 

probability  densities  were  not  explicitly  defined  in  this 

study,  they  were  described  by  determination  of  their  central 

mordents.   The  general  form  of  the  central  moment  of  a  given 

probability  density  is  given  by  equation  (11). 

R 


R 


N   (x.-p) 
E   ■ 

i  =  l 


N 


(11) 


where   pD  =  R    central  moment 
1th 


R 


1 


value  of  the  variable  X 


N   =  total  number  of  values  of  X  in  the  ensemble 

N 
M   =   E   X  .  /  N 


1=1 


1 


mean  value  of  the  ensemble 
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In  practice,  the  third  and  fourth  moments  are 
normalized  by  the  third  and  fourth  powers  of  a    (a  =   y?), 

0 

the  standard  deviation  of  the  ensemble  (Tennelces  and  Lumley, 
1972).   In  this  study  the  first  moment  is  taken  to  be  the 
mean  rather  than  zero  as  in  equation  (11). 

The  criterion  used  to  estimate  the  accuracy  of  the 
computed  data  -  the  "figure  of  merit"  -  was  the  ratio  of  a 
central  moment  of  the  computed  record  to  that  of  the  observed 
record.   This  simple  criterion  was  intended  to  show  only  two 
things:  the  percentage  of  a  given  observed  statistic  present 
in  the  corresponding  computed  record;  and  those  statistics 
which  differed  in  sign  between  computed  and  observed  records 
(resulting  in  a  negative  figure  of  merit). 

As  the    computed  data  were  compared  to  observations 

Q 

at  six   OS,  the  average  figure  of  merit  for  the  model  was 
taken  to  be  the  average  of  the  figures  of  merit  for  each  OS. 
An  average  was  found  also  for  the  three  OS  in  the  Pacific  and 
for  the  three  in  the  Atlantic, 
a.    Wind  Speed 

The  values  of  the  central  moments  and  their 
figures  of  merit  for  the  observed  and  computed  wind  speed  are 


o 

The  two  computed  records,  NCI  and  NC2,  were  compared 
to  the  observational  record  at  OSN.   In  order  not  to  over- 
emphasize these  data  when  considering  averages,  the  average 
of  the  figures  of  merit  of  NC2  and  NCI  was  used  as  the  figure 
of  merit  at  OSN  when  calculating  the  overall  average  figure 
of  merit  or  the  average  of  the  figures  of  merit  in  the  Pacific 
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tabulated  in  Table  III  and  are  plotted  against  latitude  in 
Figure  5.   The  average  computed  wind  speed  (first  moment) 
varied  between  95%  and  91%  of  the  average  wind  speed  observed 
at  the  six  OS.   Over  the  latitude  range  of  the  OS  (30°N  to 
56.5°N),  the  average  observed  wind  speed  decreased  from  a 
maximum  of  10.6  m  sec"   in  the  Westerlies  (OSP)  to  6.6  m  sec" 
in  the  North  East  Trades  (OSN)  for  the  ten  years  between  1960 
and  1969.   The  average  figure  of  merit  for  the  mean  computed 
wind  speed  is  0.92.   This  figure  is  0.91  in  the  Pacific  and 
0.93  in  the  Atlantic. 

It  is  worth  noting  again  here  that  we  are 
dealing  with  an  average  wind  speed  and  not  a  resultant  wind 
speed.   The  resultant  wind  speed  can  be  determined  by 
combining  the  average  zonal  and  meridional  components  by 
vector  addition.   The  average  wind  speed  was  determined  from 
the  wind  speed  at  each  interval  in  the  ten  year  time  series. 

The  variance  (second  moment)  of  the  wind  speed 

2    -2 
decreased  from  approximately  28  m   sec   in  the  latitudes  of 

2     -  2 

the  Westerlies  to  approximately  10  m   sec    in  the  latitudes 

of  the  Trades  -  a  latitude  dependence  similar  to  that  of  the 
average  wind  speed.   The  average  figure  of  merit  for  the 
variance  is  1.07.   Only  at  0SB  is  the  figure  of  merit  (0.89) 
less  than  unity.   The  maximum  (1.24)  is  shared  by  0SJ  and  OSN 
(NC2).   When  comparing  the  variances  of  the  computed  and 
observed  records  by  means  of  the  F  test  it  becomes  apparent 
that  any  figure  of  merit  for  the  variance  which  is  not  equal 
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to  unity  is  statistically  significant  at  the  1%  level  due  to 
the  very    large  (>6000)  number  of  degrees  of  freedom  for  each 
ensemble  (Freund,  1962).   That  is,  the  probability  that  the 
computed  and  observed  ensembles  come  from  the  same  probability 
density  is  less  than  0.01.   Relating  these  statistically 
significant  differences  in  the  observed  and  computed  records  to 
to  physical  causes  is,  for  the  most  part,  beyond  the  scope  of 
this  study.   However,  the  location  of  the  grid  point  NC2 
(at  32N)  at  the  northern  boundary  of  the  North  East  Trades 
appears  associated  with  the  high  figure  of  merit  at  this 
location.   Seckel  (o_p_  c i t . )  suggested  the  wind  speed  is  much 
more  variable  at  32N  (NC2)  than  it  is  in  the  Trades  (at  OSN), 
and  thus  a  figure  of  merit  for  wind  speed  variance  greater 
than  1.0  is  to  be  expected  between  these  records.   This  suggests 
that  the  model  can  resolve  significant  changes  in  the  varia- 
bility of  the  wind  field  which  occur  on  a  scale  of  one 
gri  d-mesh  length. 

The  standard  deviation  (square  root  of  the 
variance)  can  be  used  to  estimate  an  upper  bound  on  the 
average  fluctuation  of  the  wind  speed  from  the  ten-year  mean. 
According  to  Chebyshev's  Inequality,  the  probability  that  the 
wind  speed  will  differ  from  the  average  by  more  than  1.414 
standard  deviations  is  less  than  0.50  (Freund,  o_p_  c  i  t .  )  .   The 
value  of  1.414  a    in  the  Westerlies  was  about  7  m  sec   .   Tin's 
value  is  considerably  greater  than  the  diurnal  range  of  the 
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wind  speed  in  the  open  ocean  reported  by  Roll  (op  c  i  t .  )  . 
Thus,  most  of  the  variance  in  the  wind  speed  observed  at  OS 
over  the  ten-year  period  was  due  to  variations  of  a  period 
longer  than  twelve  hours.   The  50%  reduction  in  the  variance 
due  to  filtering  as  indicated  in  Section  II-D-2  suggests  that 
50%  of  the  variance  is  due  to  wind  speed  fluctuations  of  a 
period  greater  than  two  days. 

The  third  moment  (skewness)  of  the  computed 
wind  speed  is  higher  than  that  of  the  reported  wind  speeds 
at  all  the  OS.   The  average  figure  of  merit  for  all  the 
locations  is  1.30.   This  high  value  represents  a  tendency  of 
the  computed  wind  speed  to  exceed  reported  wind  speed 
significantly  in  cases  of  high  wind.   Of  the  locations 
studied,  this  tendency  is  most  pronounced  in  the  Pacific 
Trades.   The  fact  that  the  average  figure  of  merit  was  greater 
than  one  for  the  skewness  and  less  than  one  for  the  average 
of  the  wind  speed  suggests  that  the  overestimated  high  wind 
speeds  in  the  computed  record  were  more  than  compensated  by 
more  numerous  underestimati ons  of  low  wind  speeds.   The 
omission  of  the  effects  of  curvature,  i.e.,  the  acceleration 

of  the  wind  due  to  the  centrifugal  force,  in  the  geostrophic 

9 

equation  was  suggested  as  the  cause  of  this  effect. 


R.  Elsberry,  Dept.  of  Meteorology,  Naval  Postgraduate 
School,  personal  communication. 
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The  exaggeration  of  both  extremely  low  and 
high  wind  speeds  is  also  suggested  by  the  fourth  moment 
(kurtosis)  which  indicates  the  amount  of  extreme  data  in  the 
enseusble,  i.e.,  the  number  of  data  points  in  the  ensemble 
whicli  are  far  from  the  mean.   The  figure  of  merit  for 
kurtosis  is  highest  for  the  Pacific  in  the  transition  from 
Westerlies  to  Trades  (OSV  and  0SN-NC2).   As  can  be  seen  in 
Figure  5,  the  value  of  the  kurtosis  of  the  wind  speed 
increases  dramatically  in  the  Trades.   However,  this  increase 
is  an  artifact  of  the  normalization  by  the  fourth  power  of 
the  standard  deviation,  which  is  much  smaller  in  the  Trades 
than  in  the  Westerlies,  and  does  not  indicate  that  ^/ery    strong 
winds  occur  more  frequently  in  the  Trades  than  in  the  Westerlies 
b.    Zonal  Component 

The  first  four  moments  for  the  zonal  component 
of  the  wind  vector  are  shown  plotted  against  latitude  in 
Figure  6.   Their  values  and  figures  of  merit  are  tabulated 
in  Table  IV.   The  latitude  dependence  of  the  first  moment 
(average  or  resultant)  of  the  observed  and  computed  zonal 
components  is  similar.   The  difference  between  the  first 
moments  of  the  observed  and  computed  zonal  components  is 
greatest  at  OSP  (0.85  m  sec"  ).   While  the  average  figure  of 
merit  of  the  first  moment  is  1.0,  the  model  is  more  accurate 
in  1 1.  e  Atlantic  than  in  the  Pacific  Ocean.   In  the  Pacific  the 
figure  of  merit  ranges  from  0.43  to  1.45  and  averages  1.01. 
In  the  Atlantic  the  range  was  smaller,  from  0.89  to  1.05, 
and  the  figure  of  merit  averages  0.98. 
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This  disparity  between  Atlantic  and  Pacific 
is  evident  also  when  comparing  the  variance  (or  second  moment) 
of  the  observed  zonal  component  to  that  of  the  computed  zonal 
component.   The  figure  of  merit  ranges  between  0.83  and  0.95 
in  the  Atlantic  Ocean  and  between  0.78  and  1.02  in  the 
Pacific.   Overall,  it  averages  0.88,  indicating  that 
the  model  reproduced  about  88%  of  the  variance  of  the  zonal 
wind.   The  relatively  poor  accuracy  of  the  model  at  0SP  as 
shown  in  Figure  6  is  unexplained;  however,  the  poor  agreement 
between  0SN  and  grid  point  NC2  is  probably  due  to  the  transi- 
tion to  the  Trades  as  previously  suggested. 

The  skewness  showed  the  poorest  agreement  of 
all  the  moments  calculated  for  the  computed  and  observed 
zonal  components.   The  model  consistently  underestimated  this 
parameter;  the  figure  of  merit  averaged  only  0.63.   The  signs 
of  the  computed  and  observed  skewnesses  differed  at  0SJ ; 
however,  the  values  at  this  point  were  so  near  zero  that  the 
difference  is  insignificant. 

The  kurtosis  (fourth  moment)  of  the  computed 
zonal  component  was  consistently  higher  than  that  of  the 
observed  zonal  component.   The  average  figure  of  merit  was 
1.11,  and  was  nearly  equal  in  both  the  Atlantic  and  Pacific 
Oceans  . 

c.    Meridional  Component 

The  first  four  moments  of  the  probability 
density  of  the  computed  and  observed  meridional  components 
are  tabulated  in  Table  V  and  are  plotted  against  latitude  in 
Fi  gure  7 . 
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The  figure  of  merit  for  the  mean  of  the 
meridional  component,  averaged  for  all  locations,  was  1.22 
and  was  nearly  equal  in  the  Pacific  and  the  Atlantic  Oceans. 
The  individual  figures  of  merit  of  the  means  suggest  that  the 
computed  meridional  component  has  a  bias  of  0.3  to  0.9  m  sec" 
This  bias  is  a  fictitious  (computed)  component  from  the  south, 
which  reduced  the  average  computed  meridional  component  at 
those  stations  where  the  resultant  direction  was  from  the 
north  (OSB  and  OSN);  and  it  increased  the  meridional  component 
at  those  stations  where  the  resultant  direction  was  from  the 
south.   Such  a  consistent  bias  suggests  that  a  fictitious 
east-west  gradient  (higher  pressure  to  the  east)  in  the  sur- 
face pressure  analyses  may  occur  over  both  the  Pacific  and 
Atlantic  Oceans. 

The  variance  of  the  observed  and  computed 
meridional  components  showed  only  slightly  better  agreement 
than  was  the  case  for  the  zonal  component.   Averaged  over  the 
locations  studied,  the  model  reproduced  92%  of  the  variance 
of  the  observed  meridional  component.   Again,  OSN  (NC2) 
showed  the  highest  figure  of  merit  due  to  the  effects  of  the 
transition  to  the  Trades. 

The  skewness  of  both  the  observed  and  computed 
values  was  nearly  zero  at  all  locations  north  of  the  Trades. 
This  indicates  that,  relative  to  the  average  meridional 
component,  the  strong  positive  (northerly)  and  negative 
(southerly)  fluctuations  were  equally  distributed.   The 
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negative  value  of  the  skewness  at  OSN  indicates  that  the 
southerly  (relative  to  the  mean)  extremes  predominated. 
Because  the  skewness  was  so  near  zero  for  all  stations  north 
of  the  Trades,  figures  of  merit  at  these  locations  have  little 
meaning.   In  the  Trades,  the  model  overestimated  the  skewness 
by  up  to  98%.   This  yery    high  figure  of  merit  is  in  marked 
contrast  to  that  for  the  skewness  of  the  zonal  component 
(Table  IV).   This  discrepancy  suggests  that  in  the  Trades, 
extreme  values  in  the  calculated  record  of  meridional  com- 
ponents are  not  always  matched  in  the  observed  record. 

The  value  of  the  kurtosis  (fourth  moment)  of 
the  meridional  component  was  slightly  overestimated  by  the 
model.   The  figure  of  merit  averaged  1.19,  and  the  averages 
for  the  Atlantic  and  Pacific  Oceans  were  identical. 

The  values  of  the  skewness  and  kurtosis  of  the 
meridional  component  north  of  the  Trades,  approximately  0  and 
3.0  respectively,  suggest  that  the  values  of  the  meridional 
component  exhibit  the  well  known  Gaussian  distribution  in 
thi  s  area  . 

Of  the  four  central  moments  discussed  above, 
perhaps  the  most  important  for  oceanographi c  purposes  is  the 
variance.   This  parameter  is  a  measure  of  the  kinetic  energy 
of  the  wind,  and  is  directly  related  to  the  frictional 
coupling  of  the  ocean  and  atmosphere. 

On  the  average,  for  the  locations  studied,  the 
model  reproduced  about  90%  of  the  variance  in  each  of  the 
components  of  the  wind  vector  and  overestimated  the  variance 


of  the  wind  speed  by  7%.   In  none  of  the  cases  was  the  error 
greater  than  25%.   The  consistent  overestimati on  of  the 
variance  at  NC2  as  compared  to  NCI  indicates  that  features 
such  as  the  transition  from  the  Trades  to  Westerlies  can  be 
resolved  within  one  mesh  length  of  the  computation  grid. 

The  difference,  shown  in  equation  4,  between 
the  square  of  the  mean  wind  speed  (derived  from  individual 
pressure  analyses  as  in  this  study)  and  the  square  of  the 
resultant  wind  speed  (derived  from  monthly  mean  pressure 
analyses)  is  often  neglected  in  studies  of  large  scale  air- 
sea  interaction  (Namias,  op  c  j t .  ) .   As  can  be  seen  from 
Tables  III  through  V,  this  difference  as  observed  at  OSP  is 

The  resultant  wind  squared  is  only 


92.5  m   sec" 


20.7  m   sec"  .   In  this  case  the  resultant  square  is  less 
than  20%  of  the  mean  squared  (112.2  m2  sec   ). 

The  difference  between  the  ten-year  mean  and 
ten-year  resultant  is  of  course  larger  than  the  corresponding 
difference  over  a  month  due  to  the  longer  period  (seasonal 
and  year-to-year)  variation.   However,  this  example  does  show 
that  the  difference  between  the  square  of  the  mean  and  the 
square  of  the  resultant  can  not  be  considered  insignificant 
when  dealing  with  the  energetics  of  large  scale  air-sea 
interaction. 

2  .    Ten-Year  Trend 

The  e  s  t  i  m at  ion  of  a  time  series  of  wind  conditions 
using  pressure  analyses  from  different  sources  raises  a 
problem  cited  by  Bakun  ( op  c  i  t .  )  .   That  is,  any  long-term 
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trend  in  the  wind  record  may  be  due  to  different  methods  of 
analyzing  the  pressure  distribution  in  successive  time 
intervals.   For  example,  if  improved  analysis  techniques  or 
more  numerous  pressure  reports  allowed  better  resolution  of 
tight  pressure  gradients  there  could  be  a  fictitious  increase 
in  the  wind  speed  derived  from  the  improved  pressure  analyses 
Obviously,  such  a  fictitious  increase  should  not  be  reflected 
in  the  long-term  trend  of  observed  wind  speeds. 

In  order  to  investigate  this  possibility,  the 
ten-year  changes  in  the  observed  and  computed  wind  speeds  and 
components  were  calculated  using  equation  (12)  (Bendat  and 
Persol,  1971). 


u(t)  =  b_  +  b,t 


10 


where  u(t)  =  value  of  trend  line  at  time  t,    and 


(12) 


N 


N 


2  ( 2 N  + 1  )   E  un    -  6   E  nu, 
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h  N  (N-l )  ( N  + 1 ) 


(14) 


The  sample  interval,  h,  in  equation  (14)  is  assumed 
to  be  equal  to  1.   However,  the  gaps  in  the  record  were 
not  filled  prior  to  the  determination  of  the  ten-year  trend 
Because  the  vast  majority  of  the  gaps  consisted  of  a  single 
missing  data  point,  this  should  have  little  impact  in  the 
determination  of  the  ten-year  change. 
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value  of  the  ensemble 
N  =  total  number  of  values  in  the  ensemble. 

The  change  over  the  ten-year  period,  a,  is  given  by  equation 

(15): 

A   =  bjN  .  (15) 

The  values  of  A  for  the  wind  speed  and  each  of  the  components 
are  shown  plotted  against  latitude  in  Figure  8. 

In  the  Westerlies  the  reported  wind  speed  has 
decreased  over  the  ten  years  between  1960  and  1970.   In  the 
Trades  (OSN)  it  has  increased  "jery    slightly.   The  decrease 
appears  to  be  most  pronounced  in  the  Atlantic  Ocean  with  a 
maximum  decrease  of  1.16  m  sec"   at  OSJ.   For  the  other 
locations  the  ten- year  change  was  less  than  0.6  m  sec"  . 
Computed  wind  speeds  showed  changes  similar  to  those  observed 
at  all  locations  with  the  exception  of  OSP,  where  the  slight 
increase  of  the  computed  speed  (in  opposition  to  the  observed 
decrease)  is  unexplained. 

With  the  exception  of  OSJ  and  OSN  the  values  of  the 
observed  zonal  components  show  a  general  increase  over  the 
ten  years  of  the  study.   It  should  be  noted  that,  due  to  the 
negative  sign  of  the  prevailing  zonal  component  in  the 
Westerlies  and  the  positive  sign  in  the  Trades,  this  repre- 
sents a  decrease  in  the  prevailing  zonal  components  at  all 
the  locations  except  at  OSJ. 
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At  OSJ  the  observed  increase  in  the  prevailing 
zonal  component  was  not  accompanied  by  an  increase  in  the 
prevailing  meridional  component.   In  fact,  the  prevailing 
observed  meridional  components  decreased  at  all  locations, 
with  the  exception  of  OSD. 

Although  there  are  noticeable  differences  between 
the  trends  of  the  computed  and  observed  series,  i.e.,  OSV 
(zonal),  OSD  (meridional)  and  OSP  (speed),  the  similarity  in 
latitude  dependence  and  magnitude  between  observed  and 
computed  records  suggests  that  the  trends  in  the  computed 
series  are  largely  controlled  by  processes  controlling  the 
trend  of  the  observed  series,  and  are  not  due  to  artificial 
changes  in  the  pressure  data. 

3 .    Correlation  and  Regression  Analysis 

This  section  deals  with  the  linear  correlation  and 
non-linear  regression  analyses  carried  out  on  the  ensembles 
of  paired  (computed  and  observed)  data.   Non-linear 
(parabolic)  regression  was  done  in  order  to  determine  the 
non-linearity  of  the  relation  between  the  computed  wind  data 
and  those  observed  at  the  OS. 

The  linear  correlation  coefficient,  R ,  is  found 
using  equation  (16): 


N 

E  ( X  -  x  )  ( y  -  y  ) 


N        9  N        , 

i     (*n-*>   *  (yn--y) 

n=l         n=l    n 


(16) 


40 


where  x  =  computed  value 
n 

y  =  observed  value 


N   xn 

*   (4) 
n  =  l 


N 


N   yn 

s  (if) 

n=l 


N 


The  value  of  R   is  a  measure  of  the  dependence  of  the  varia- 
bility of  one  variable  on  the  variability  of  the  other  through 

linear  coupling  (Panofsky  and  Brier,  op  c  i  t .  )  .   Table  VI 

2 
indicates  the  values  of  R  and  R   for  the  wind  speed  and 

components,  and  the  number  of  paired  values  upon  which  the 

correlation  is  based.   In  all  cases,  the  correlation  between 

computed  and  observed  values  of  the  components  exceeds  that 

between  computed  and  observed  speed.   Averaged  over  the  six 

OS,  the  model  explained  65%  of  the  variance  in  each  of  the 

components  and  42%  of  the  variance  in  the  wind  speed. 

The  lower  correlation  for  wind  speed  can  be 

attributed  to  the  combined  effect  of  errors  in  both  the  zonal 

and  meridional  components.   If  the  computed  wind  vector  is 

W  ,  then 
c  • 


W   =  Z   i  +  M  j 
c    c      c 


(17) 


where  Z   and  M   are  the  computed  zonal  and  meridional  com- 

c      c  ' 

ponents  and  i,  j  are  unit  vectors  in  zonal  and  meridional 

directions.   If  we  also  define  the  wind  error  vector  W   as: 

e 
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A  A 


We  -    Hc  -  W0  , 


(18) 


where  W   is  the  observed  wind  vector 
o 


Wo  =  V  +  V  ' 


(19) 


we  have,  in  terms  of  components, 


W  =  ( Z   -  Z  )  i  +  ( M   -  M  )  j  ; 
e     c    o        c    o 


(20) 


and  the  wind  speed  error  (magnitude  of  the  error  vector)  is 
given  by 


W   =  {(Z   -  Z  )2  +  (M   -  M  )2}2 
e     v  c    oy    v  c    o' 


(21) 


Thus,  the  error  (and  resulting  lack  of  correlation)  in  the 
wind  speeds  should  always  be  greater  than  or  equal  to  the 
error  (and  lack  of  correlation)  in  the  components;  i.e., 
the  correlation  coefficient  for  the  wind  speed  should  always 
be  less  than  that  for  either  of  the  components. 

On  the  average,  the  computed  data  correlated  best 
with  observations  for  the  zonal  components.   In  the  Atlantic 
Ocean  67%  of  the  observed  variance  in  the  zonal  component  was 
explained  by  the  model.   However,  the  meridional  component 
showed  higher  correlation  than  the  zonal  component  for  five 
of  the  seven  series. 
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In  addition  to  linear  correlation,  non-linear 
regression  between  the  computed  and  observed  records  were 
also  determined.   The  regression  equations  are  given  below: 


Y  =  A   +  A,X  +  A0X' 
O     1       2 


(22) 


X  =  B   +  B,Y  +  B9Y' 
o     I     d. 


(23) 


where   Y  =  computed  value 
X  =  observed  val  ue 

The  coefficients  were  found  by  solving  the  follow- 
ing set  of  simultaneous  equations  (24a  -  24c)  (Spiegel,  1961) 
for  the  computed  data  as  a  function  of  the  observations,  and 
solving  the  inverse  set  (replacing  Y  by  X  and  solving  for 

B  ,  B ,  and  B0)  for  the  observations  as  a  function  of  the 
Oi       2 


computed  data 


ZY    =  A  N  +  A, EX  +  A0EX' 
o      1       2 


EXY   -  A  EX  +  A,EX2  +  A0EX3 
o       I        2 


IX2Y  =  A  EX2  +  AnEX3  +  A0ZX4 
o        i        2 


(24a) 
(24b) 
(24c) 


Table  VII  lists  the  coefficients  A  ,  A, ,  and  A0  for  equation 

0    1        2       ^ 

(22).   Table  VIII  gives  the  coefficients  for  the  inverse 
equation  (23).   The  values  of  N ,  the  number  of  paired  data 
values  in  each  series,  are  found  in  Table  VI. 
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Comparing  A,  to  Ap  and  B,  to  B2  indicates  that  an 
essentially  linear  relation  exists  between  the  computed 
records  and  observed  records.   The  non-linear  coefficients, 
A2  and  B2,  are  two  orders  of  magnitude  less  than  A,  and  B, 
for  all  series  except  for  the  meridional  component  at  OSN. 
At  this  location  the  non-linear  term  was  only  one  order  of 
magnitude  less  than  the  linear  term. 

Figures  9  and  10  are  plots  of  the  two  regression 
equations  for  wind  speed  at  OSJ  and  OSN  (NC2).   These  loca- 
tions showed  the  highest  and  lowest  correlation  between 
observed  and  computed  wind  speed.   The  separation  of  the  two 
regression  curves  is  due  to  the  scatter  of  the  computed  and 
observed  wind  speed  data  being  compared,  as  shown  in  Figures 
11  and  12  (together  with  Table  IX).   The  scatter  of  the  data 
is  at  least  partially  due  to  neglecting  some  of  the  signifi- 
cant physics  in  the  model,  i.e.,  stability,  and  the  curvature 
of  the  isobars  in  the  pressure  analyses.   The  area  where  the 
regression  lines  converge  (the  area  of  least  scatter)  is  the 
area  of  best  fit  of  the  model.   For  OSJ  this  region  is 
centered  around  11  m  sec"  ;  for  OSN  it  centers  around 
7  m  sec-  .   In  all  cases  this  region  of  best  fit  is  near  the 
mean  value  of  the  series. 

If  the  relation  of  the  computed  parameter  to  the 
observed  parameter  changes  with  time,  then  no  single  time- 
independent  equation  can  express  this  relation.   For  the 
purposes  of  describing  the  relative  accuracy  of  the  model  at 
various  locations,  however,  the  relation  between  the 


44 


observed  and  computed  data  is  assumed  to  be  shown  by  a  single 
"true"  regression  line.   The  "true"  regression  line  is  defined 
as  that  line  bisecting  the  distance  between  the  two  regression 
curves  X  =  f(Y)  and  Y  =  f(X).   This  distance  is  measured  per- 
pendicular to  the  X=Y  line.   The  "true"  regression  lines  were 
found  by  graphical  means  for  the  wind  speed  at  OSJ  and  OSN 
and  are  plotted  in  Figures  9  and  10. 

These  figures  indicate  that  a  consistent  feature 
of  the  model,  overestimati on  of  high  observed  wind  speeds 
and  underestimation  of  low  values,  is  most  pronounced  in 
the  low  latitudes  (OSN)  where  observed  wind  speeds  below 
about  11  m  sec"   are  underestimated.   This  underestimation 
occurs  to  a  lesser  degree  at  OSJ  up  to  wind  speeds  of  about 
14  m  sec"  .   The  scatter  cannot  be  reduced  by  applying 
correction  factors  (deduced  from  the  regression  analysis)  to 
the  output  of  the  model;  however,  taking  average  values  of 
the  model  output  over  periods  longer  than  12  hours  will 
provide  a  more  accurate  time  series  (at  the  expense  of 
reduced  time  resolution).   Th.is  is  a  result  of  the  elimination 
of  the  extreme  values  in  the  time  series  by  the  averaging 
process.   As  the  averaging  interval  increases,  the  mean  over 
that  interval  approaches  the  mean  of  the  time  series  in 
question.   This  set  of  means  will  then  exhibit  less  scatter 
than  the  original  series  as  they  will  be  in  the  neighborhood 
surrounding  the  point  where  the  two  regression  lines 
[X=f(Y)  and  Y=f(X)]  intersect  the  X=Y  line. 
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4.    Error  Analysis 

Two  sources  of  error,  in  addition  to  the  neglected 
physics,  are  separation  distance  and  the  effects  of  friction 
in  low  latitudes.   Investigation  of  the  accuracy  of  the 
surface  pressure  distribution  from  which  the  wind  records 
were  calculated,  a  third  possible  error  source,  is  beyond  the 
scope  of  this  thesis. 

a.    Separation  Distance 

The  computational  grid  points  do  not  exactly 
coincide  with  the  locations  of  the  OS  where  the  observations 
were  made.   The  distance  between  the  grid  point  and  OS  is 
defined  as  the  separation  distance.   The  effect  of  separation 
distance  on  the  mean  error,  i.e.,  the  difference  between  the 
mean  of  the  observed  record  and  that  of  the  computed  record, 
is  shown  plotted  against  separation  distance  in  Figure  13. 

The  mean  error  of  the  wind  speed  shows  no 
relation  to  separation  distance.   However,  such  a  relation 
is  recognizable  for  the  components  of  the  wind  vector.   The 
meridional  component  mean  error  shows  a  slight  increase  with 
increasing  separation  distance.   Although  there  is  consider- 
able scatter,  there  does  not  seem  to  be  any  significant 
difference  between  Atlantic  and  Pacific  OS  in  this  regard. 

This  is  obviously  not  the  case  for  the  zonal 
component  however.   The  Atlantic  OS  show  little  mean  error 
in  the  zonal  component.   At  Pacific  OS,  however,  a  rather 
pronounced  dependence  on  separation  distance  is  evident. 
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This  dependence  suggests  the  following: 

(1)  The  ten-year  average  of  the  north-south 
pressure  gradient  (from  which  the  zonal  wind  is  derived)  is 
more  accurately  described  in  the  North  Atlantic  than  in  the 
North  Pacific  Ocean. 

(2)  The  accuracy  of  this  average  gradient 
within  200  km  of  OS  in  the  Pacific  is,  in  part,  dependent  on 
the  distance  from  the  OS. 

In  addition  to  the  effect  of  separation 
distance  on  mean  error,  its  effect  on  the  correlation  between 
computed  and  observed  records  was  also  investigated.   The 
results,  shown  in  Figure  14,  indicate  no  correlation  coeffi- 
cient dependence  on  separation  distance  comparable  to  that 
for  mean  error. 

b.    Friction  in  Low  Latitudes 

Low  latitude  wind  fields  play  a  very  important 
part  in  large  scale  air-sea  interaction  (Berjknes,  1969; 
White  and  Walker,  1973).   However,  the  wind  fields  in  this 
region  are  typically  difficult  to  model  because  they  are 
strongly  ageostrophic.   As  no  data  comparable  to  OS  observa- 
tions were  available  from  low  latitudes  for  direct  verifica- 
tion, this  section  is  included  as  an  attempt  to  evaluate  the 
accuracy  of  the  model  in  this  region. 
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In  a  recent  study,  Brummer,  ejt^  aj_.  (1974) 
summarized  data  suggesting  that  the  strong  ageostrophic 
component  of  the  surface  wind  vector  in  low  latitudes  was  due 
to  the  effects  of  friction.   These  effects  resulted  in  a 
reduced  ratio,  C,  of  surface  wind  speed  to  geostrophic  wind 
speed  and  an  increase  in  the  deflection  angle,  a.      Table  X 
lists  the  values  of  C  and  a  as  reported  by  Briimmer,  e_t  al  .  , 
and  those  used  in  the  model  at  corresponding  latitudes. 
(The  values  of  a    as  used  in  the  model  are  for  surface  wind 
speeds  of  2.5  m  sec"  .)   Although  the  values  of  C,  as  tabu- 
lated by  Brummer,  e_t  aj_.  ,  show  no  latitude  dependence  similar 
to  those  used  in  the  model,  they  agree  fairly  well  in  magni- 
tude.  The  value  of  a  used  in  the  model  were  close  to  those 
tabulated  by  Briimmer  ,  et  aj_.  ,  with  the  exception  of  a  at  2°N. 

At  this  latitude  the  values  of  a    and  C  as  used 
in  the  model  could  result  in  an  error  in  the  orientation  of 
the  surface  wind  vector,  and  corresponding  errors  in  the 
components.   For  example,  in  the  case  of  pure  zonal  flow,  an 
error  of  30°  in  a  would  reduce  the  zonal  component  to  0.866 
of  its  original  value  and  introduce  a  fictitious  meridional 
component  equal  to  0.5  of  the  original  zonal  value.   Such 
a  fictitious  meridional  component  in  the  computed  record 
could  lead  to  erroneous  conclusions  concerning  the  correlation 
of  zonal  and  meridional  wind  stress  anomalies.   However,  for 
the  latitudes  north  of  about  5°N,  the  parameterization  of  a 
and  C  appear  to  be  reasonable. 
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B.    SPECTRAL  ANALYSIS 

The  method  of  analysis  described  in  the  preceding 
section  allowed  comparison  of  computed  and  observed  values 
only  as  ensembles.   This  section  describes  the  results  of  a 
comparison  of  frequency  components  constituting  the  observed 
and  computed  time  series  using  spectral  analysis.   No  attempt 
will  be  made  to  discuss  the  atmospheric  dynamics  implied  by 
the  spectra. 

1 .    Autocorrel ati  on 

The  autocorrelation  function  is  the  basis  of 
spectral  analysis.   It  describes  the  relative  amplitude  of 
the  different  frequency  components  present  in  a  given  time 
series  (Bendat,  Peirsol,  op  c i t  ) .  The  autocorrelation  Rvv(x) 

X  A 

can  be  defined  as: 


xx 


(t)  =  lira  j     j     {X(t)x(t+t)>  dt 


(25) 


where 

t  =  lag  i  nterval 

T  =  total  time  of  the  record 

t  =  time  variable 

For  this  study  T  consisted  of  the  ten  years  from 
1960  through  1969,  corresponding  to  7  306  twelve  hourly  incre- 
ments of  t.   The  lag  interval  ranged  between  0  and  one  year 
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(731  twelve-hour  increments).   Because  of  the  gaps  in  the 
observed  record  and  the  filtering  performed  on  the  time 
series  (see  Section  1 1  - D )  ,  the  spectral  analysis  data  are 
presented  only  for  frequencies  greater  than  0.5  cycles  day"'. 

Table  XI  presents,  at  selected  frequencies,  the 
values  of  the  autocorrelation  function  of  the  computed  and 
observed  time  series  at  OSB  and  OSV.   The  frequency  is  deter- 
mined by  the  lag  interval  x,  which  is  the  period  of  the 
frequency  component  in  question.   The  discussion  of  the 
autocorrelation  function  is  intended  to  illustrate  the  nature 
of  the  periodicity  of  the  observed  time  series  and  the 
differences  between  them.   (The  accuracy  of  the  model  at 
specific  frequencies  is  best  described  by  the  phase  and 
coherence  functions  discussed  in  Section  II-B-3.)   All  of  the 
computed  series  were  expected  to  show  autocorrelation  higher 
than  the  observed  data  for  short  lags  (0  to  ^8  hours)  due 
to  the  averaging  nature  of  the  geostrophic  approximation 
(Murk,  1960).   This,  however,  was  not  the  case  at  either  OSB 
or   OSV.   However,  because  the  data  were  heavily  smoothed  at 
these  frequencies  the  differences  may  be  masked.   The  values 

of  R   (t)  for  the  observed  and  computed  wind  speed  at  OSB  and 
xx  r 

OSV  ere  shown  plotted  against  lag  time  in  Figures  15  and  16. 
The  annual  variation  in  wind  speed  is  more  pronounced  at  OSV 
than  at  OSB  in  spite  of  the  fact  thai.  OSB  is  at  a  higher 
latitude.   The  reduced  values  of  R   ( t )  at  annual  and  semi- 

A  A 

annLil  frequencies  at  OSB  can  be  attributed  to  a  stronger 
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random  component  in  the  wind  speed  at  these  frequencies. 

(OSV  perhaps  is  influenced  by  the  Monsoon  regime  as  suggested 

by  Malkus,  1962.  ) 

The  nature  of  Ryy(t)  for  the  zonal  component  differs 

J\  A 

dramatically  between  OSB  and  OSV  as  shown  by  Figures  17  and 
18.   At  OSB  the  zonal  component  is  essentially  random  for 
lags  longer  than  4  days.   This  is  not  the  case  at  OSV  where 
annual  variation  is  pronounced.   The  flattened  minimum  at 
OSV  (Figure  18)  indicates  that  the  semi-annual  periodicity 
(between  100  and  250  days)  of  the  zonal  component  was  reduced 
at  this  location  by  random  fluctuations. 

The  meridional  component  at  both  OSB  and  OSV, 
Figures  19  and  20,  show  predominately  random  nature  for  lags 
greater  than  4.0  days.   However,  Figures  19  and  20  show  "chat 
the  meridional  component  at  OSB  has  a  slightly  more  periodic 
nature  than  at  OSV. 

The  computed  series  at  both  OSB  and  OSV  show  the 
same  general  features  as  do  the  corresponding  observed 
records.   In  most  cases  even  the  minor  details  of  R  v(t)  for 

A  A 

the  observed  record  are  also  present  in  the  R   (x)  for 

X  X 

computed  record.   This  suggests  that  the  aperiodic  as  well  as 
the  periodic  nature  of  the  wind  is  accurately  reproduced  by 
the  model. 

2 .    Energy  Density 

For  a  given  frequency,  f,  the  energy  density, 
G   (f)  can  be  expressed  in  terms  of  the  autocorrelation 

A  A 

function  as  foil ows  : 
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Gxx<f>  ■  2   /  Rxx^  e_i^fT  ** 


(26) 


where  RYY(x)  =  autocorrelation  for  lag  t.   This  relation  is 

A  A 

a  Fourier  transform  from  period  space  (x)  to  frequency  space 
(f).   The  area  under  a  given  segment  of  the  energy  density 
curve  is  proportional  to  the  variance  of  the  time  series  in 
that  frequency  range. 

The  energy  density  at  frequencies  equivalent  to  the 
lags  (periods)  listed  in  Table  XI  in  the  preceding  section  are 
listed  in  Table  XII.   A  general  underestimation  of  the  energy 
by  the  model  is  indicated  by  this  table  and  by  the  reduced 
value  of  the  variance  of  the  computed  component  records  as 
indicated  in  Section  ( 1 1 1  - A ) .   Most  of  the  "missing"  variance 
in  the  computed  records  occurs  at  frequencies  less  than 
0.1432  cycles  day"   (periods  greater   than  1  week).   The 
dependence  of  the  energy  density  on  frequency  for  the  computed 
and  observed  records  is  shown  in  the  upper  diagrams  of  Figures 
21  through  26.   It  should  be  noted  that  the  plot  of  the  loga- 
rithm of  the  energy  density  tends  to  overenhance  the  absolute 
differences  at  high  frequencies  and  to  mask  them  at  low 
frequencies . 


The  scale  of  the  spectral  energy-density  plot  varies 
from  figure  to  figure  as  a  result  of  attempting  to  maximize 
the  resolution  of  data,  which  varied  significantly  from  plot 
to  plot. 
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However,  the  plots  indicate  that  the  model's 
accuracy  is  not  significantly  dependent  on  frequency  and  that 
the  geographical  differences  of  the  observed  wind  spectra 
are  accurately  reproduced  in  the  computed  wind  spectra. 
Geographical  differences  are  especially  evident  in  the  spectra 
of  the  components;  this  indicates  that  for  purposes  of  large 
scale  numerical  modeling  of  sea-air  interaction,  uniform, 
zonal -averaged  wind  fields  cannot  be  considered  realistic. 
3.    Coherence  and  Phase 

Coherence  is  a  measure  of  the  correlation  between 
two  time  series  at  a  given  frequency;  as  such,  it  is  the  best 
estimator  of  the  accuracy  of  the  model.   It  can  be  defined 
in  terms  of  the  energy-density  spectra  of  each  of  the  two-time 
series  of  two  variables,  x  and  y,  and  their  cross  spectra. 
The  equation  for  coherence,  y x v s  is  given  in  equation  (27). 


xy 


{Gxv(f)r 
(f)  ■       y 


Sxx(f)  Gyy(f) 


(27) 


where  G   (f)  and  G   (f)  are  the  energy  densities  of  the 

computed  and  observed  time  series  as  described  in  the  previous 

section.   The  cross-spectral  density  function,  G   (f),  is 

xy 

defined  in  terms  of  the  cross  correlation  function  between 
the  variables  x  and  y,  R   (  t  )  (Bendat  and  Peirsol,  o_p_  ci  t .  )  . 


G   (f)  =  2 

xy v  ' 


j       R_(Oe-i2ltfT  dx 


xy 


(28) 
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where  the  cross  correlation  function  is  given  by 


xy 


(t)  =  Tim   J  x(t)y(t+x)dt 


(29) 


"!"-><» 


If  v   (f)  =  1,  the  records  are  perfectly  correlated  at 
xy 

frequency  f;  for  Yyv(f)  =  0>  no  relation  exists  between  x 

xy 

and  y  at  that  frequency. 

Another  estimator  of  the  accuracy  of  the  model 
derived  from  the  cross-spectral  density  function  is  the 
cross-spectral  phase  angle,  0   (f).   This  function  measures, 
at  frequency  f,  the  average  angular  difference  between  the 
time  series  of  x  and  y,  and  is  defined  in  terms  of  the  real 
and  imaginary  parts  of  the  cross -spectral  energy -den si ty 
function  : 


0   (f)  =  tan"1  [Q   (f)  /  C   (f)] 
*xyx  '  xy *  '    '       xy v  ; J 


(30) 


where  Q   (f),  the  imaginary  part  of  G   (f),  is  called  the 

quadrature  spectral-density  function;  and  C   (f),  the  real 

part  of  G   (f),  is  called  the  coincident  spectral- density 

function.   C   (f)  and  0   (f)  can  be  defined  as  follows: 
xy v  '  xy 


C   (f )  -  h    (G   (f )  +  G   (f ) } 
xy  xy       y  x 


(31) 


Q   (f )  =  \    {G   (f )  -  G   (f )} 
wxyv/    2xy^;     yxv/ 


(32) 
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For  this  study  a  positive  phase  angle  indicates 
the  computed  record  lags  the  observed  record  in  time.   For 
example,  a  cross-spectral  phase  angle  equal  to  +9°  at  a 
frequency  of  0.1  cycles  day"   indicates  that  the  component 
of  computed  record  at  this  frequency  lags  the  observed  record 
by  0.25  days  or  6  hours. 

The  values  of  the  coherence  and  phase  (cross- 
spectral  phase  angle)  at  OSB  and  OSV  for  the  frequencies  of 
interest  are  tabulated  in  Tables  XIII  and  XIV.   In  Figures 
21  through  26  the  coherence  and  phase  are  plotted  in  the 
lower  two  diagrams.   At  both  OSB  and  OSV  there  is  a  gradual 
increase  in  coherence  with  decreasing  frequency.   As  suggested 
in  Section  III-A-3,  the  coherence  between  computed  and 
observed  wind  speed  is  lower  than  is  the  coherence  between 
the  observed  and  computed  components.   At  OSB  there  is  a 
constant  increase  in  0xv(f)  with  increasing  frequency.   For 
the  wind  speed  and  the  components  of  the  wind  vector  this 
phase  difference  is  approximately  45°  at  f=0.5  cycles  day"  , 
indicating  the  computed  record  lags  the  observed  by  about  6 
hours . 

For  OSV  there  is  no  consistent  phase  difference 
between  the  observed  and  computed  records  similar  to  that  at 
OSB;  however,  the  coherence  is  similar  to  that  for  OSB,  and 
approaches  unity  at  periods  of  six  months  and  a  year. 
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The  data  presented  in  this  section  suggest  that 
this  model  reproduces  the  major  features  of  the  observed 
energy  density  spectra  of  wind  speed  and  orthogonal  components 
for  frequencies  between  0.5  and  0.0027  cycles  day"   (periods 
of  12  hours  to  one  year).   It  also  suggests  that  spectra  for 
the  components  of  the  wind  vector  cannot  be  deduced  from  the 
corresponding  spectra  of  the  wind  speed. 


IV.   SUMMARY 


A.    ACCURACY  OF  THE  MODEL 

The  results  of  this  study  suggest  the  following: 

1.  The  relation  between  the  computed  and  observed 
records  of  wind  speed  and  components  of  the  wind  vector  is 
essentially  linear. 

2.  The  model  exaggerates  both  high  and  low  values  of 
wind  speed  and  components  and  this  exaggeration  is  most 
pronounced  at  low  latitudes. 

3.  The  accuracy  of  the  model  as  shown  by  correlation 
with  observations  is  slightly  higher  in  the  Atlantic  compared 
to  the  Pacific  Ocean;  and  it  is  higher  in  mid-latitudes  than 
in  either  high  or  low  latitudes. 

4.  The  ten-year  (1960-1969)  trend  of  wind  conditions 
between  29°N  and  56.5°N  is  similar  for  both  the  computed 
records  and  the  observed  records;  and  observations  indicate 
a  slight  reduction  (-0.5  m  sec"  )  in  wind  speed  over  the 
northern  hemisphere  between  35 N  and  5 ON. 

5.  The  model  explains  between  30%  and  70%  of  the 
variance  of  the  parameters  measured  at  12- hour  intervals. 

6.  There  is  a  fictitious  southerly  component  incorpo- 
rated in  the  computed  records  at  all  grid  points,  possibly 
the  result  of  a  fictitious  east- west  gradient  in  the  pressure 
analyses  used. 
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7.  The  deflection  angle  between  the  surface  wind  and 
the  geostrophic  at  latitudes  between  0°  and  5°N  as  used  in 
the  model  is  substantially  less  than  the  values  thus  far 
observed. 

8.  The  accuracy  of  mean  wind  conditions  derived  from 
the  model  increases  with  the  length  of  the  averaging  interval, 
as  indicated  by  increasing  coherence  with  longer  period 
fluctuations. 

B.    POSSIBLE  APPLICATIONS 

The  last  conclusion  points  to  a  possible  application  of 
the  model  in  addition  to  its  current  operational  use.   Monthly 
means  of  wind  conditions  derived  from  12-hourly  surface 
pressure  analysis  using  the  model  should  exhibit  even  greater 
accuracy  than  do  the  12-hourly  values.   The  mean  wind  para- 
meters derived  in  this  way  would  be  true  averages  and  not 
resultant  wind  parameters.   A  historical  record  of  this 
nature  should  permit  a  more  realistic  quantitative  approach 
to  sea-air  interaction  studies  than  has  been  the  case  to  date. 
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TABLE  I 


Monthly  distribution  of  surface  pressure  analyses.   The 
number  of  analyses  used  is  in  the  numerator,  and  the 
denominator  equals  the  number  of  twelve  hour  intervals  in 
the  month. 
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TABLE  IX 


Scatter  diagram  symbols  (for  Figures  11  and  12) 
and  the  corresponding  number  of  coincident 
computed  and  observed  wind  speed  pairs. 
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TABLE  X 

Ratio  of  surface  wind  speed  ( V s )  to  geostrophic 
wind  spaed  ( V g )  and  deflection  angle  (in  degrees) 
at  low  latitudes  as  tabulated  by  Brummer,  ejt  al  . 
and  as  formulated  in  the  model. 
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TABLE  XIII 


Coherence  between  observed  and  computed  time  series 
for  Ocean  Station  Bravo  and  Ocean  Station  Victor 
at  selected  frequencies. 


PERIOD 
(DAYS) 

FREQ. 
(DAI  ~1) 

BRAVO 

VICTOR 

ZONAL 

meridJ 

SPEED 

ZONAL 

1-JERiD. 

SPEED 

2.0 

.5000 

.705 

.501 

.562 

.430 

.686 

.624 

4.0 

.2500 

.790 

.750 

.692 

f780 

.856 

.780 

7.0 

.1432 

.855 

*806 

.760 

.742 

.814 

.522 

29.6 

.0338 

.927 

.812 

.705 

.943 

e<*888 

.783 

185.0 

.0054 

.943 

.927 

o8S9 

i  .970 

.804 

.731 

370.0 

.0027 

.920 

.946 

.963 

I  .978 

.770  | 

.937 

wmmmmmmmtmam 

TABLE  XIV 


Phase  difference  (in  degrees)  between  observed  and 
computed  time  series  for  Ocean  Station  Bravo  end 
Ocean  Station  Victor  at  selected  frequencies. 


PERIOD 
(DAYS)  . .. 

FREQ. 
(DAY  ~1) 

BRAVO 

VICTOR 

ZONAL 

MERID. 

I 
SPEED 

hZSML. 

MERID „  I  SPEED 

2.0 

.5000 

0.3 

29  o  2 

39*9 

-9.8 

-8^2        2.6 

4.0 

.2500 

8C0 

29.0 

14.9 

-7.2 

4,5 

5.6 

7.0 

.1432 

14.2 

7.1 

-12.3 

-19.7 

-4.9 

-4.4 

29.6 

.0338 

1.7 

0.1 

11.2 

0.9 

12.0 

-1.5 

185.0 

.0054 

-0.6 

-3i2. 

~206 

-3.7 

-3.4 

-0.2 

370.0 

.0027 


-5.6 

-5.9 

-3.5 

-1.7 

-15.1 

8.2 
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Figure  2.  Latitude  dependence  of  the  deflection  angle 
for  various  wind  speeds  (W  in  m  sec-'')  as  formulated 
in  the  model . 
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Figure  3.   Latitude  dependence  of  the  deflection  angle 
as  found  L»y  Roll  (broken  line),  Carstensen  (solid 
line),  and  as  formulated  in  the  model  for  wind  speeds 
( W )  of  5  m  sec-1  and  15  m  sec-1  (dashed  lines). 
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Figure  4 .   Result  of  five  point  binomial  filter  applied 
to  a  typical  portion  of  a  zonal  component  time  series 
Original  series  (solid  line)  and  filtered  series 
(dashed  line)  arr>  shown  in  the  upper  graph.   The 
residual  series  is  shown  in  the  lower  graph. 
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Figure  8 .   Total  change  (integrated  linear  trend) 
over  the  ten  year  period  from  January  1960  thro 
December  1969  in  the  observed  (solid  line)  and 
computed  (dashed  line)  wind  speed  and  zonal  and 
meridional  components  (in  m  see-!)- 
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Figure  11.   Scatter  of  observed  and  computed  wind 
speeds  (in  n  s e c - 1 )  at  Ocean  Station  Juliet. 
Tabic  IX  defines  the  symbols. 
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Figure    12.       Scatter    of    observed    and    computed    wind 
speeds    (in    m    sec"^)    at    Ocean    Station   Hovember. 
Table    IX    defines    the    symbols. 
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Figure  13.   Mean  error  (in  m  sec"  )  dependence  on 
separation  distance  for  wind  speed  arid  the  zonal 
and  meridional  components.   Estimated  best-fit 
lines  are  shown  for  components. 
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Figure  14.   Linear  correlation  coefficient  (R) 
dependence  on  separation  distance  for  wind  speed 
and  the  zonal  and  meridional  components. 
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Figure    15.       Autocorrelation    of    observed    (CO)    end 
computed    (XX)    wind    speeds    at    Ocean    Station    bravo 
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Figure    16.       Autocorrelation    of    observe-*!    (CO)    and 
computed    (XX)    wind    speeds    at    Ocean    Station    Victor 
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Figure    17.      Autocorrelation    of    observed    (CO)    and 
computed    (XX)    zonal    components    at   Ocean    Station 
Bravo . 
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Figure    18.       Autocorrelation    of    observed    (CO)    and 
computed    (XX)    zonal    components    at    Ocean    Station 
Victor. 
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Figure    19.       Autocorrelation    of    observed    (CO)    and 
computed    (XX)    meridional    components    at    Ocean 
Station    Bravo. 
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Figure    20.      Autocorrelation    of    observed    (CO)    and 
computed    (XX)    meridional    components    at    Ocean 
Station    Victor. 
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Figure  21.   Log  energy  density  (in  in   sec"  ),  coherence,  and 

phase  difference  (in  degrees)  for  computed  and  observed 
wind  speeds  at  Ocean  Station  Bravo. 
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Figure  24.   Log  energy  density  (in  m   sec   ),  coherence,  and 

phase  difference  for  computed  and    observed  zonal  components 
at  Ocean  Station  Victor. 


'■  . 


> 

I— 

I — I 
LO 

z:    •". 

!_U     <? 

a 

> 

CC 
SjJ 
ZI     m 

LU      4 


d 

CC 

^— 
01    ^, 


0.300 


U-  4l,w 


0.500 


z: 

CC 
LU 
X 

D 
U      o 


'^^Y^^^ 


I 
0.100 


-1     7  -.  -. 


300 


0.500 


Frequency  (cycles  /  day) 


LO 

en 


000        '".':;.         0.200        0.J00        0-400        0. 


O.'IOO 


2     -  2 
Figure  25.   Lon  energy  density  (in  m   sec   ),  coherence, 

and  phase  difference  for  computed  and  observed  meridional 

components  at  Ocean  Station  Bravo. 
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Figure  26.   Log  energy  density  in  m   sec   ),  coherence  and 
phase  differences  for  computed  and  observed  meridional 
components  at  Ocean  Station  Victor. 
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